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'p^i Abstract. We discuss harmonic generation in the case of laser field-dressed 

Hamiltonians that are invariant under so-called dynamical symmetry operations. 
Examples for such systems are molecules which exhibit a discrete rotational symmetry 
of order A'^ (e.g. benzene with A^ = 6) interacting with a circularly polarized laser 
^ I field and single atoms in a bichromatic field, with the two lasers having circular 

O ' polarizations. Within a general group theory approach we study the harmonics one 

c/3 , obtains from the interaction of a laser pulse and a circular molecule. When the system 

^^' is in a pure field-dressed state the known selection rule fcA^ ± 1, A; = 1, 2, 3, . . . results. 

f^ . However, other lines are observed when recombinations with states of a symmetry 

different from the initial one become important. This is the case for realistic laser pulses 
(i.e., with a finite duration). In particular when the fundamental laser frequency (or one 
^ , of its multiples) is resonant with a transition between field-dressed states. Numerical 

'^r ' ab initio simulations, confirming our analytical calculations and illustrating the power 

^ii^ ■ of the group theory approach, are presented. 
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1. Introduction 

The generation of harmonics by atoms and molecules in a laser field is a topic that 
has been thoroughly studied during the last fifteen years (see [Q for recent reviews). 
The big interest lies in the potential of atoms and molecules in laser fields as sources 
of short-wavelength and short pulse radiation. It is well known that when a strong 
linearly polarized laser field interacts with an atom only odd harmonics are generated. 
In more complex systems a detailed study of the symmetry transformations is required 
for predicting which harmonics are allowed or not 0, ^. In a recent work 0] it was 
shown that peculiar selection rules can be obtained for molecules which are invariant 
under a discrete rotational symmetry Cat (like benzene where A^ = 6). When this 
kind of molecules interacts with a circularly polarized laser pulse of infinite duration, 
the selection rules for the harmonics are n = kN ± 1, where k G A/^. and n is the 
harmonic order. This class of molecules is particularly interesting because they can 
work as a filter, i.e., the first allowed harmonic is n = A^ — 1, and therefore, if A^ is large, 
the emitted frequency can be very high. All systems whose Hamiltonian is invariant 
under a "dynamical rotation", i.e., a discrete spatial rotation plus a time translation, 
generate selection rules of the kind kN ± 1. Together with cyclic molecules another 
example of such a system is the single atom interacting with two-color laser fields with 
circular polarizations and frequencies one an integer multiple of the other. In both the 
molecular and the atomic case it is possible to obtain these selection rules also from the 
conservation of the angular momentum component perpendicular to the polarization 
plane. 

In a previous paper [Q we have presented results from numerical simulations for 
a benzene model-molecule showing that, when a realistic, finite pulse is taken into 
account, other lines than those expected from the kN ± 1 selection rule can be present. 
In this paper we will show that, if embedded in a general group theory treatment, 
also this case can be understood in a systematic way. Our analysis will allow us to 
achieve a full understanding of the molecule harmonic spectrum. Moreover, we will 
show that the field-dressed level scheme can be inferred from the harmonic spectrum. 
The particular features and differences of a finite pulse compared to an infinite one will 
be also addressed. All the laser fields considered are circularly polarized in the xy-plane, 
and atomic units are used throughout the paper. We apply the dipole approximation 
which is excellent for the laser frequencies and intensities under consideration. 

2. Dynamical Symmetries 

The Hamiltonian H^o of a cyclic molecule interacting with a single laser field and the 
Hamiltonian i^at of a single atom interacting with two circularly polarized lasers are, 
from a symmetry point of view, equivalent. Both are invariant under certain dynamical 
rotations. We will consider only a single electron active. This assumption does not 
introduce any loss of generality to our derivations |0]. Therefore H^o in cylindrical 
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coordinates (p, ip, z) reads 

+ ^P cos(^-^t), (1) 

where £ is the amphtude of the electric field. If V^mo(p, V', z) is a model potential for a 
cyclic molecule with A^ ions the Hamiltonian H^o is invariant under the transformation 

i 

A A 27r 27r \ 

In the case of the single atom in the two-color field if at is 

I d ( d\ 1 92 -92 

£ £ 

+ ^P cos{ip-ujt) + ^p cos{ip + r]ujt), (3) 

where £1 and £2 are the electric field amplitudes of the two laser fields and opposite 
polarization has been assumed. The second laser has a frequency that is rj times the 
frequency of the first laser, with 77 integer. It is nice to see that it is really the second 
laser which generates the particular symmetry. In fact, all the terms except the last one 
of if at, are invariant under the continuous transformation 

Poc= U^^ + 0,t^t+-\, (4) 

where 6 has, so far, an arbitrary value. If no other term is present in ifat no harmonics 
are generated. Forcing also the last term to be invariant under such a transformation 
leads to 

271 A- 
e + ri9 = {ri + l)e = 2nk^e = . (5) 

Therefore the Hamiltonian iiat is invariant under the discrete transformation 

and this is equivalent to the case of a cyclic molecule with A^ = 77 + 1. Hence the selection 
rules for this case are equal to those derived in for the cyclic molecules, n = kN ± 1. 
The case where the two lasers have the same polarization can be handled with the same 
method, leading to similar selection rules. Further details will be discussed elsewhere 
IP. In the following the two lasers are always considered with opposite polarization. 

An experiment with the two frequencies u and 20; was performed [^, and its 
theoretical aspects were recently discussed within a semiclassical three-step model 0. 
Selection rules for the case of an atom interacting with a circularly polarized field 
of frequency u and a linearly polarized field of frequency A^o; were also derived |§], 
confirming numerical results obtained previously B. 
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In the derivation of the kN ± 1-selection rule one assumes that only a single non 
degenerate, bound state is involved in the harmonic generation process. With this 
restriction the selection rules can be also derived simply with a change of reference 
frame. In the rotating laser field frame the observer "sees" an ionic potential oscillating 
with a frequency Nu in a static field; in the two-color atomic case instead, after the 
change of frame, the atom is at rest in the static field and it experiences an oscillating 
field of frequency Nuj with N = rj + 1. Hence, in this frame the dipole has only Fourier 
components of frequency kNu with k G A/". Going back to the laboratory frame where 
the molecule is at rest, the Fourier components are shifted ±uj and the selection rules 
presented in P] are obtained. Keeping the assumption of a single bound state involved 
in the process, it is, moreover, possible to derive the same selection rules also arguing 
about the angular momentum only. Let us first discuss about the cyclic molecule: 
because of the symmetry of the molecule, an expansion of the ground state in spherical 
harmonics Yf^ has to look like ^o = YlT=o 'l2i>kN '^k^i^t^ + c.c, i.e, the ground state is 
expanded over functions that have the required azimuthal symmetry. The wavefunction 
is then invariant under a rotation of 271 /N . Let us assume that the laser field has a 
polarization a^ , i.e., the projection of the angular momentum of a photon on the field 
propagation direction is equal to 1. If a harmonic of polarization a^ is emitted, the 
molecule, for keeping the symmetry, hence has to absorb kN + 1 photons so that the 
final change in the projection of the angular momentum is kN . On the contrary, if the 
harmonic emitted has polarization cr~, the molecule must absorb kN — 1 photons. The 
consequence is that for each k the higher harmonic kN + 1 has the same polarization 
as the laser field, and the lower kN — 1 the opposite one, as it was obtained with the 
dynamical symmetry method in . For the atom the derivation is similar [0, [Tl| . Let 
us assume an electron is in an s state (the quantum numbers / and m are zero). Let the 
first laser of frequency u have a polarization a^ and the second laser of frequency rj uj 
a polarization a". If a harmonic a^ is emitted the sum of all the components along z 
of the angular momentum carried by absorbed photons has to be +1. If from the first 
laser the atom absorbs p photons, it must absorb p — 1 from the second laser. The total 
absorbed energy is poj + {p — l)rjuj = uj[p{l + rj) — rj] = uj[{p — 1)(1 + rj) + 1]. Since 
p is arbitrary, with k = p + 1 and A^ = r^ + 1 we see that the {kN + 1)*^ harmonic 
is emitted. With the same argument, starting from a harmonic of polarization a~ , we 
obtain kN — 1. 

3. Group Theory 

From a group theory point of view the Hamiltonian H^^ of a cyclic molecule without 
the laser field has a symmetry described by the group Vi\fh |]12|, ^ . When the molecule 
interacts with a laser field which has no temporal envelope its total Hamiltonian, given 
by (p, can be written as ifmo = -f^mo + W{t) where W{t) oc pcos{ip — ut) is the 
molecule-laser interaction. 

First we study the symmetry group of the total Hamiltonian -ffmo- H^^ has the 
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symmetry Vjyh and W is described by Qoo, i.e, an extension to infinity of tlie so called 
"cyclic dynamical symmetry group" 0, defined by Ql = {Pl, Ph ■ ■ ■ , Pl~^ , Pt — 1} 
where L is the order of the group and the number of its elements. Pl is the operator 
already introduced in (^ which acts in space and time, and is defined as 

^^^+_i^« + _j. (7) 

Combining V^h and Qoo we have that the symmetry group of the Hamiltonian H^^ is Qn 
where A^ is the number of ions in the molecule which determines the discrete rotational 
symmetry. 

Because Qn is a cyclic group, it is also Abelian, and therefore it has the following 
properties: (i) each element forms a class by its own and hence there are A^ irreducible 
representations, (ii) each irreducible representation is one-dimensional, i.e., the elements 
are complex numbers. The group Qn is isomorphic to the well known group Cn and 
therefore the same character table apply. It is interesting to look how, from a physical 
point of view, the symmetry of the molecular Hamiltonian changes when the laser field 
is present. The group V^h contains AN elements, i.e., N rotations about an axis of the 
iVth order, N rotations through an angle vr about horizontal axes, N reflections a^ at 
the vertical planes and N rotary-reflection transformations C%(Jh [ll2[- Looking at the 



field-dressed Hamiltonian ([1|) it appears that only the rotations about the vertical axis 
are still symmetry transformations. In fact, the reflection at the vertical planes and the 
rotations about horizontal axes would require a space-dependent time transformation. 
The reflection ct/j at the horizontal plane would change the polarization of the laser field. 
Hence from the 4A^ possible operations we had for H^^ only A^, i.e., the rotations about 
the vertical axis, keep the Hamiltonian if mo invariant. 

For a certain group Cn the type of its irreducible representations depends on A^. 
If A^ is odd there are one real and A^ — 1 complex representations, if A^ is even there 
are two real and N — 2 complex representations. Half of the complex representations 
is the complex conjugate of the other half. Usually in the literature, the two complex 
conjugated representations are considered equivalent and indicated with E, normally 



reserved for two-dimensional representations [|T2|. In our case, because the Hamiltonian 
(HI) is not time- reversal invariant, these representations are not equivalent. Therefore 
we have to keep all the A^ representations separated. 

We call each of the A^ irreducible representations R^ where i is an index between 
and A^ — 1. The generator of the i-th. representation is given by e'^'^^/^. Therefore, it 
follows that the (m + l)-th entry of the i-th representation is e^'^'^^"^/^ (the first entry 
is given by the identity). Because of the unidimensionality of the representations the 
product of two representations, Rq, and R/3, is the irreducible representation Ra+i3- If 
a + P > N, Ra+f3 is equivalent to Ra+i3-N- Thus each entry of the product table is again 
a representation R^ with < £ < A^ — 1. 

Three conclusions can be drawn now: (i) from the definition of the representation 
generator we have that the complex conjugated representation of R^ is Rat-^, (ii) the 
representation Rq has all characters equal to 1 (this is the well known "totally symmetric 
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irreducible representation") and it is present for every A^, (iii) when A^ is even RAr/2 is 
the other real representation. 

The description given so far in terms of group theory will prove particularly useful 
later on when we analyze the transitions between the field-dressed states of the molecule. 
Indeed, whether an electronic transition between two states, each of them corresponding 
to an irreducible representation, is allowed or not can be determined examining the 
product of the symmetries of the initial state, the inducing operator, and the complex 
conjugated of the final state. If the product leads to the totally symmetric irreducible 
representation the transition is possible, i.e., Rgn (^ Rop (S) ^mi = Ro in our notation. 

4. Harmonic generation 

Provided the laser pulse is sufficiently adiabatic and the emitting molecules can be 
considered uncorrelated it is a good approximation to calculate the harmonic spectrum 
from the Fourier-transformed dipole expectation value. With the wavefunction \E'(r, t) 
we can write the probability to get the nth harmonic as 

e'""*dt /"^*(f,t)pe^^'^*(f,t)df 



x{n) oc n 



n^ 




^*(f, t)pe^("'^*^'^)^(f, t)df dt 



(8) 



For a system whose Hamiltonian is periodic in time with a period r = 27r/u; we can 
apply the Floquet theorem (see, e.g., \T^\)- The wavefunction of a certain state of the 
system can be written as \E'j(r, t) = e~^^'*$j(r, t) where C,i is the so-called quasi-energy 
and ^i{f,t + T) = $,j(r, t). The functions $j(r, t) are eigenfunction oiHmo = [Hmo — ^-§i], 
i.e.. 

In what follows we refer to $i(r, t) as a Floquet state. If the system is in a pure Floquet 
state we can derive the selection rules as presented in 0. However, due to pulse shape 
effects or other perturbations the system might be not in a pure Floquet state, but 
rather in a superposition of them, 

^{f,t) = J2(3^e-'^''Mr,t). (10) 

i 

Here, for simplicity, we assume the /3jS constant (or slowly varying) in time. The quasi- 
energies are given by ^i = Aj — iF, where the real part Aj is the field-dressed energy of 
the state and the imaginary part Fj corresponds to half the ionization rate. The latter 
would lead to finite harmonic line widths of Lorentzian shape. In what follows we will 
neglect this effect keeping only the real part of C,i for our symmetry analysis. Eq. (|^) 
then becomes 



x{n) oc n'^ 



j2 lAr / / $.pe^("-*^^)<i>, dfdt 
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i>j 



We introduce A^ = Aj — Aj and the following operator 



:ii) 



G 



(iJ) 



^i({nuj-Al)tT'p) 



Noting that the functions ^i{f,t) form an extended Hilbert space |T^ Eq. (^) can be 
rewritten as 



x(n) oc n'^ 



j2m'm\G^''^m) + j2p:PAmGi''¥j)) 



i>j 



(12) 



where the double brackets stand for integration over time and space. 

We can think about the different operators Gzp as operators which induce transitions 
where the initial and the final states can be the same or different. Such transitions are 
accompanied by the emission of harmonics. If the system is stable in one Floquet state 
the expansion in ([T^) reduces to only one term. In particular, when we consider the 
system being in the Floquet ground state ^o{f,t), x{^) is proportional to the Fourier- 
transform of the expectation value of the dipole calculated with the Floquet ground 
state, i.e., 

2 



{i,i)l 



xin)ocn^ ((<l>o|G^'^'|$o)) 



(13) 



From (p!2|), the structure of the dipole spectrum can be derived. Because Pj^^Pj^^ = 1, 

{($,|G^'^')|<I>,)) = {{PM<I>^\PNG'^''^PN'\PN<^J)), (14) 

and writing the operator G^'"' explictly we have 

,i{(ni^-Ai)tT'p) 



pNG^i'^pN^ = Pjvpe'«"-"^^)*^^)P^1 



pe 



eit("-i^Ti) 



G 



(»j) 



.iTf("-l^Tl) 



(i,i) 



Applying M times P/v to the operator GV we obtain 



M r<{i,j) \ f)-liM 



[p^r G^'^^[p, 



N 



G 



{i,j) 



.iT^("-^Tl) 



(15) 



(16) 



From a group theory point of view the term in square brackets in (|T^) is simply the 
generator of the corresponding irreducible representation. From (|T^ we see how the 
harmonic generation operator behaves under the dynamical rotations of the symmetry 
group of the Hamiltonian. Given the irreducible representations of the states $j(r, t) 
and ^j{r, t) we can derive to which representation the operator G-^' must correspond in 
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order to satisfy Eq.(|T^). In representation notation, we must have R* (^ Rop (^ Rj = Rq, 



where Rop is the irreducible representation of the operator G-^'"' . From the definition 
of the character table of an Abelian group |jl2[ we know that Pn^j (^, t) is equal to the 
phase e^~^v~ times ^j{f,t) and similar for ^i{f,t). Also PnG^ P^ generates a phase. 
For allowing the emission of a harmonic the product of these three phases must be the 
real unity. Consequently, the operator G^ has to belong to the representation R^^-^j. 
Imposing the term in square brackets of Eq. ([15|) to be equal to the generator of the 
representation R^^-^, leads to 



; 2tt I 



'^i^-,^ .27r(fi-<?,0 



e^wi--=^Ti) ^ gi^^:^^ ^ ^(^ _ M ^ 1) = '^^(^^^-^j) + 2vrfc. (17) 
From this follows 

G^'^^ — ^ n = kN + ii + ^±l (18) 

where £*■ = ii — ij. Introducing the extended index £j = £j — ^, i.e., the symmetry 
index minus the field-dressed energy of the state in frequency units, we can rewrite (^) 
in the final form 

G^'^'^ — > n = kN + ii±l (19) 

where £*• = ii — ij. 

Varying n we force the generic operators G=p to behave under the dynamical 
rotations according to the different representations R^ of the symmetry group, and for 
each i certain harmonics are emitted. For the special case G^' one obtains n = kN ± 1, 
reproducing the results in [0]. They have been calculated taking £* = and A^ = 
because the initial and the final state are the same. G^' behaves according to the 
totally symmetric representation Rq. 

It is worth to be noted that if the energy difference of the field-dressed states 
corresponding to the irreducible representations R^- and R^^ is equal to uj{ii — if) we 
have a kind of accidental degeneracy: the lines corresponding to such transition cannot 
be distinguished from those arising from the single state processes. 

The operators G=p can behave according to any representation. Hence all couplings 
between the Floquet states of the molecule are in principle possible. This of course does 
not mean that they have equal probability to occur. Considering that the ground state 
corresponds to Rq, and the other R^ represent excited states we can derive the following 
symmetry properties for the relevant single transitions: (i) Rq — > Rq: symmetry Rq, [ii) 
Ro -^ Re'- symmetry R^, (Hi) R^ — > Rq: symmetry R^ = Rn-i = R-e, {iv) Re — ^ R^: 
symmetry Rq. 

It is evident that when the frequency uj (or one of its integer multiples) of a finite 
laser pulse is nearly resonant with a transition to a certain state, this state becomes 
more easily accessible, and therefore transitions involving that state are more likely to 
happen. On the other hand, when the laser frequency is tuned off from all the resonances 
mostly the field-dressed ground state, i.e. Rq, is involved and it will be the only one 
giving a significant contribution to harmonic generation. 
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The fact that selection rules for different initial and final states contain the terms 
Aj is particularly interesting for the case where the rotating field removes a degeneracy. 
The measurement of the distance between the satellite lines arising from the two states 
gives an estimation of the field-dressed level-splitting. 

5. Benzene Molecule 

The benzene molecule has a symmetry V^h, but as long as we consider only functions 
which are all even or odd with respect to reflection at the molecule plane {ah in the 
standard notation) it is sufficient to consider only the symmetry operations of the 
reduced group Vq. The character table of the group Vq is given in Tab. 1. 

Connecting the four lowest states of the benzene molecule and the irreducible 
representations of the group Vq, i.e., Ai, Bi, E2, Ei leads to the following 
energetic sequence: Ai (ground state), Ei (first excited, two-dimensional, i.e., two- 
fold degenerated), E2 (second excited, two- fold degenerated), Bi (third excited, non- 
degenerated) |jl3[- In the ground state Ai and Ei contain two and four vr-electrons. 



respectively, and the other levels are empty ||T6[. As it was discussed in the previous 
section, when the circularly polarized laser field has an infinite duration, the total 
Hamiltonian H^^ = H^^ + W would belong to the dynamical symmetry group Qq 
which is isomorphic to Cg. The character table of the groups Cq and Qq is shown in 
Tab 2. For the group Cq the standard notation is reported, for the group Qe instead we 
use the notation that has been previously introduced. With the laser field switched on 
the degeneracies connected to the representations Ei and E2 are removed. If we take 
for example the degenerated level Ei of the field-free system, in an infinite pulse it is 
reducible, and knowing the irreducible representations of the group Cq we can easily 
decompose it into irreducible parts: Ei — > E^ + E\ = Ri -|- R5. The same argument 
can be applied to the level E2, E2 — > Eg + Eg = R2 + R4. 

What is worth to be stressed is that, contrary to the standard case in quantum 
mechanics textbooks, here the two representations which are complex conjugated, like 
Ri, R5 and R2, R4, correspond to two different energy levels and therefore the so-called 
"physically irreducible" representations really generate two independent states. This 
particular feature is due to the non-invariance of the total Hamiltonian H under a time 
inversion. A similar case is obtained when a magnetic field is considered. 

At this point, with an oscillating field infinite in time, we can think about the 
field-dressed molecule states in pure terms of Rj, where i runs from to 5. This system 
is stable, i.e., there are no transitions between the states corresponding to different 
representations, unless external perturbations are present. The selection rules in this 
case can be derived from the product table. Tab. 3. When the pulse has a finite length 
we can consider the pulse shape as a perturbation of a system that has a Vq symmetry. 
A transition between two states is possible only if the operator inducing the transition 
has the "proper" symmetry. As it can be derived from Tab. 1 the representations Ai, 
El, E2 and Bi correspond to rotations of angle 0, |, ^ and vr, respectively. In terms of 
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£ they correspond to i = 0, £ = 1, i = 2 and £ = 3. 

6. Numerical Simulations 

In this Section we illustrate and test our group theoretical method by applying it to 
harmonic spectra obtained with the help of ab initio numerical simulations. In order to 
keep the numerical effort feasible we restrict ourselves to a two-dimensional (2D) model 
where the molecule plane and the rotating electric field are properly represented. We 
study a single active electron in a ring-shaped potential with Cn symmetry (A^ = 6 in 
the benzene case). Results from a fully 3D calculation could be quantitatively different 
but the structure of the spectrum would be the same. We use the potential [Q 

V{p, if) = -—==£=== [a cos (iV ^) + 2 - a] (20) 

where ro is the radius of the molecule (ro = 2.6), and p and ip are polar coordinates. /3 is 
the parameter which gives the "degree of smoothness" of the potential and determines 
the width of the ground state along the ring, a moves the local maxima of the potential 
keeping the minima constant (the parameter a is introduced in order to avoid the 
presence of points where V^ = for finite p, because that could generate non-physical 
ionization). Finally, A is the "strength" of the potential. For our simulations we chose 
a = 0.99 and (3 = 0.38. We used A for tuning the ionization energy of the molecule. 
For more details about the model and its properties see |Q. 

When the correct ionization energy for the benzene molecule is taken, the four 
states are relatively weakly bound and hence, as it has been already discussed in [Q, for 
inducing a multiphoton process with a low ionization rate a low frequency is required. 
As the main aim of the numerical simulations presented in this work is to test ab initio 
the analytical derivations discussed in the previous sections, we chose A = 1.0 which is 
less demanding for what concerns the run-time because a higher frequency can be used 

i- 

If with R we indicate the field-free energy of the level corresponding to the R 
representation of the T>q group, with A = 1.0 the energies of the four states are: 
Ai = -1.292, El = -1.226, Ea = -1.098 and Bi = -1.027, and the energy gaps 
between the excited states and the ground state are Qi = 0.0658, ^^2 = 0.1937 and 
Q3 = 0.2648. In order to have a convenient way of referencing we introduce the notation 



L!j! for indicating a transition between the state Rj and Kj and its corresponding line 



in the spectrum. L is the £ which gives the transition symmetry. The lower index 
distinguishes between the two signs in equation (|T8|). The upper indices specify the 
initial and the final state, respectively. The order of the states in the upper index gives 
information about the sign of the line shift due to the distance between the states. 
Most of the non-single-state transitions discussed here connect an excited state with 
the ground state. Transitions between excited states are possible but much less likely 
to occur. The position of the lines due to a "transition" from the ground to an excited 
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state is L^! = kN + i — Aj ± 1, and for the opposite process L!^ = kN — i + Aj ± 1, 
where Aj is the distance between the ground state and the state Rj. 

When the laser is switched on the states move because of the dynamical Stark effect 
and therefore Aj is given by Ai = Qi + Si, where fii is the separation between the two 
states when no field is present and 6i is the relative shift due to the presence of the 
oscillating field. 

At this point, as we know the selection rules and the values of the field-free energy 
intervals Qi, we can predict the structure of a spectrum where also non-single state 
processes are involved. The only uncertainty is due to the values of the field shifts 
6i. However, those shifts can be determined by performing several simulations with 
different laser intensities: when the field is low 6i is small (in the limit of zero field 
Aj — > Qi) and the lines are located in the proximity of the position expected from the 
unperturbed level scheme. With increasing field the shift of the states can be followed. 
Those shifts can become relatively large but the interpretation of the spectrum remains 
always unambiguous. 

An example of a spectrum, where also lines occurring from processes between 
different states are present, is shown in Fig. 1. There, a sine-square pulse of 40 cycles 
and frequency uj = 0.0942 interacting with the model molecule has been simulated. The 
peak field strength S was 0.14. It is possible to figure out for each peak the transition 
which generates it. The same symbol indicates the same process but with different 
multiplier k, i.e., different number of photons. A certain process is repeated every 
N = 6 frequency units. The states that together with the ground state play a role are 
the first and the second excited, i.e., Ei, and E2. The first state generates transitions of 
type 1, the second state of type 2. With the parameters used for the simulation of Fig. 
1 the role of the second state is more important than that one of the first excited state. 
This is mainly related to how strong the resonance with a particular state is. With 
a different frequency the relative intensity of the lines connected to the two processes 
can differ strongly. From the fact that the degeneracies of the states Ei and E2 are 
removed one could expect more lines than present in Fig. 1. Let us discuss for example 
about the states E2, that in our case are the most important role players. Each of 
the states E2 for a given k generates two lines given by 2^ = kN -|- 2 — A2 ± 1 and 
two given by 2^° = kN — 2 -|- A2 ± 1. Therefore eight lines are expected. However, 
using E2 and E2 to distinguish the two states E2, we find the shifts A2 = fig -|- ^2 and 
A2 = ^2 + ^2- Considering that ^2 — ^2 i^ follows that if the difference of the shifts 
6i due to the laser field is very small the eight lines are gathered in four couples, each 
of them containing two very close lines, difficult to resolve. Moreover, what makes the 
distinction between the two lines in a couple possible is not the absolute value of 6i but 
the difference between the two shifts. As the two states without field are degenerate, 
this difference is the splitting generated by the field. With our notation we have for the 
first couple of lines, e.g., 2_ — 2_ = ^2 ~ ^2- The same applies for the other three 
couples. This explains why in Fig. 1 only four lines due to the states E2 are observed. 
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In order to obtain a larger splitting, more easily detectable, we have to study 
how the shift 6i behaves with respect to the different field parameters. It was found 
that with lower frequencies the separation among the states which without field would 
be degenerate becomes more significant, as one would expect from a ponderomotive 
scaling. Moreover, it has been observed that the value of each Si can become relatively 
large compared to fii. In Fig. 2 we show how the location of certain lines moves when 
the field is increased. 

The dipole spectrum obtained with the lower frequency u = 0.0785 (all other 
parameters the same as before) is shown in Fig. 3. Each of the lines in Fig. 1 is splitted 
in two, confirming our predictions. 

In order to investigate the role of pulse shape effects simulations using trapezoidal 
pulses with different rampings have been performed. In Fig. 4 two spectra for pulses 
of the same length but with all cycle and 3 cycle ramping, respectively, are reported. 
It appears that in both case the overall structure is like in the case of the sine-square 
pulse. It is worth to note that when the ramping is shorter. Fig. 4b, the relevance of the 
"extra" lines with respect to the O-lines originating from the single state transitions is 
higher than in the case of longer ramping. Fig 4a. This is due to the fact that when the 
ramping is shorter the pulse is less adiabatic and therefore transitions between different 
Floquet states are induced more likely. In the limit of an infinitely long ramping only 
the lines would be present. 

7. Conclusions 

We presented a general group theoretical approach to harmonic generation by systems 
possessing a discrete rotational symmetry Cn, in particular a ring-shaped molecule in a 
circularly polarized laser field (e.g., benzene). In the simple case where the initial and 
the final field-dressed electronic states are the same the known selection rule kN ± 1, 
A; = 1, 2, 3, . . . was obtained. We demonstrated that this selection rule follows also from 
a change of reference frame as well as angular momentum conservation in a straight 
forward manner. Since due to this selection rule less harmonics are allowed within a 
fixed frequency range such systems might be beneficial for generating short wavelength 
radiation more efficiently, i.e., without wasting a vast amount of laser energy in undesired 
laser harmonics. However, we showed that in general other harmonics different from the 
expected ones at the positions {kN±l)uj are present as well. Those lines carry important 
information about the laser field-dressed level scheme of the molecule since energy 
differences between different states are involved. Transitions between field-dressed states 
occur because the laser pulse has a temporal shape instead of being infinite as assumed 
in the derivation of the kN ± 1-selection rule. For the case of the ring-shaped molecule 
we showed that our group theoretical method is capable to predict also the position of 
those extra lines. This was demonstrated with the help of harmonic spectra obtained 
by an ab initio numerical simulation of a benzene model molecule interacting with a 
finite laser pulse. The extra harmonics are particularly pronounced in the cases where 
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the laser frequency (or a multiple of it) becomes resonant with electronic transitions 
between field-dressed states. Harmonic peak splittings and shifts were fully understood 
in the framework of the group theory approach. From the symmetry point of view the 
case of a cyclic molecule in a circularly polarized laser field is equivalent to a single atom 
in a two-color field with one frequency an integer multiple of the other. Of course, the 
groups representing the unperturbed systems are different. 
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Table 1. Character table of the group V^ 
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Table 2. Character table of the group Qq and its isomorphic group Cg. 
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Table 3. Product table of the irreducible representations of the group Cg when the 
degeneracy due to the time inversion is removed. 
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Figure 1. A part of the dipole spectrum for a benzene model molecule. The laser 
frequency was 0.0942. Together with the lines expected by the single state approach 
(X and Y) other lines due to recombinations with states of different symmetries are 
present. The second excited state plays an important role. 
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Figure 2. Dipolc spectra for different laser peak intensities: £ = 0.10 (solid line), 
£ = 0.08 (dotted), and £ = 0.14 (dashed). The laser frequency was 0.0942. The shift 
of the lines due to the ac Stark effect can be relatively strong. The lines move right or 
left according to which was the initial state. 
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Figure 3. Dipolc spectrum of the benzene model molecule at a lower laser frequency 
Lu — 0.0785 where the peak splitting due to the removal of degeneracies is more clearly 
visible. 
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Figure 4. A part of the dipole spectra for two trapezoidal pulses of equal total length, 
42 cycles, but different ramping. In (a) the ramping was 11 cycles and in (b) 3 cycles. 
As expected, when the ramping is shorter, the lines due to transitions between different 
states are enhanced with respect to those involving a single state, 0_ and 0+. 



